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Abstract 

We discuss the diffuse emission of gamma rays and neutrinos from galaxy clusters 
in the viable models for structure formation in the universe. We use a self-consistent 
picture for cluster formation and evolution starting from a primordial density per- 
turbation spectrum, and a realistic modelling for the distribution of the intergalac- 
tic medium which is abundantly present within galaxy clusters. We find that an 
evolving population of clusters can produce a fraction ~ 0.5 -j- 2% of the diffuse 
gamma-ray background (DGRB) observed by EGRET. This result is robust and 
is weakly dependent on the cosmological scenario and on the degree of evolution 
of the inter galactic medium (IGM) in distant clusters, because the bulk of the 
sources contributing to the DGRB is located at redshifts z ^> 0.2. We also found a 
correlation between the non-thermal, gamma-ray and the thermal X-ray emissions 
from these structures. Using this result, we derived a list of gamma-ray clusters 
observable with the next generation 7-ray detectors. Finally, we briefly discuss the 
possible relevance of galaxy clusters for neutrino astronomy and for very high energy 
particle astronomy. 



1 Introduction 



The EGRET experiment on board CGRO revealed the existence of a dif- 
fuse gamma-ray background (hereafter DGRB) at the level of Idgrb = 9-6 • 
10 _7 .E£^ 1±ao5 cm^s-hr^GeV- 1 [1] in the energy range 0.03 + 10 GeV. 
However, a recent reanalysis of the EGRET data [2] found that the level of 
the DGRB is systematically lower by a factor ~ 20% in the energy range 
~ 0.1 -f- 4 GeV. The DGRB is observed at high galactic latitudes b > 10 deg 
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and such an evidence suggested an extragalactic origin for this diffuse back- 
ground. Nonetheless, the specific origin of the DGRB is still under debate. In 
fact, the EGRET experiment [3] has a poor angular resolution (9 min ~ 1 deg) 
so that it is hard to discriminate among different origins of this extragalactic 
background. Specifically, it is still difficult to discriminate between a purely 
diffuse nature of the DGRB (see e.g. [4]) and the option of a DGRB made by 
a superposition of unresolved, discrete sources. 

The large number of identified AGNs and flat spectrum radio quasars (here- 
after FSRQ) in the EGRET sky ([5] [6] [7]) suggested that most of the DGRB 
can be produced by a non-resolved population of AGNs, the actual fraction 
of the DGRB produced by FSRQ and BLLacs being in the range ~ 40 -j- 95% 
(see e.g. [8]). Separately, it has been evaluated that a fraction ~ 42 -h 97% 
of the DGRB could be ascribed to blazars [9]. However, the flatness of the 
spectrum of the DGRB seems to favour the possibility that BLLacs could be 
the major contributors to the DGRB of extragalactic origin [10]. Erlykin et al. 
[73] reviewed the various AGN contributions and quoted that the fraction of 
the DGRB produced by the observed AGNs is ~ 65%. The DGRB fractions 
previously reported may be subject to a revision (~ 25% increase) if the recent 
reanalysis [2] of the EGRET data is adopted. 

On account of the large theoretical uncertainties and of the present observa- 
tional precision of the EGRET detectors, it is still hard to discriminate among 
the different proposed possibilities, even though a fluctuation analysis of the 
EGRET data should give more precise indications on the nature and origin of 
the DGRB. 

Beside the discrete, unresolved source case pictured for the origin of the 
DGRB, there have been some pioneering works [11] [12] [13] [14] suggesting 
that a relevant fraction of the DGRB could be produced by extended sources 
through hadronic collisions of cosmic ray (hereafter CR) protons interacting 
with the protons of the Inter Galactic Medium (hereafter IGM) which is abun- 
dantly present within galaxy clusters (see [15] for a review). In this alternative 
picture, the CR's are assumed to be produced within clusters (we will discuss 
in Sect. 4 some of the possible sources) where also a population of protons and 
electrons is residing in the form of a hot (with temperatures T ~ 10 7 -t- 10 8 K), 
tenuous (with electron number densities n e ~ 10~ 3 cm -3 ), chemically enriched 
and massive (with mass fractions Mjgm/M ~ 0.05 -r 0.3) plasma: the IGM. 

The proposed mechanism has an essential ingredient in the confinement of the 
CR's within clusters where they are produced; this point, already realized by 
some authors [13] [14], is responsible for the net increase in the probability of 
interaction per proton with respect to the case of a straight line propagation. 
The increase factor can be estimated to be ~ ct c i/R c i ^ 600, where t c i ^ Hq 1 is 
the age of the cluster and R c i is its size. Cosmic rays produced within a cluster 
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during all its lifetime can thus produce gamma rays through the production 
and the subsequent decay of neutral pions: 

P + p^7l° + X, 7T°^7 + 7. (1) 

Note that in the same interactions, charged pions are also produced, which 
determine a neutrino emission through the following channels: 

p + p-t^ + X, 71-* -> t^VfSpn), -> e ± + P M (i/ /i ) + u e (P e ) .(2) 

We will also discuss the relevance of these last processes in Section 7 below. 

Using the gamma ray production from clusters of galaxies according to eq. 
(1), Houston et al. [11] suggested that the total extragalactic gamma ray 
intensity detected above 35 MeV [16], J 7 5.5 ■ 1CT 5 cm~ 2 s~ l sr~ l , could 
be ascribed, for a large fraction, to galaxy clusters. They predicted a level 
Zy ~ 5 • 10 _5 CTO _2 s _1 sr _1 at energies above 35 MeV, assuming an observed 
local cluster space density, n c \ ~ 7.3- 10~ 5 Mpc~ 3 , integrated out to the Hubble 
radius, Rh = 6-10 3 Mpc, and neglecting any cosmological effect. More recently, 
Dar & Shaviv (hereafter DS [12]) reanalyzed the problem in the light of the 
EGRET data [1] and calculated the contribution to the DGRB from CR inter- 
actions in the intracluster gas, under the assumption that the energy density 
of CR's in clusters is the same as in our own galaxy (universality). With this 
assumption, Dar & Shaviv [12] predicted a level J 7 (> 100 MeV) « 1.2 • 10~ 5 
photons cm" 2 s" 1 sr _1 , which could explain the whole amount of the DGRB of 
extragalactic origin. In a following paper, Berezinsky, Blasi & Ptuskin (here- 
after BBP [13]) relaxed the ad hoc assumption of universality, and estimated 
the CR energy density in clusters due to various possible sources of CR, using 
the condition of diffusive confinement of CRs. In their approach, BBP [13] 
showed that it is impossible to fulfill the universality condition with the usual 
CR sources in clusters, emphasizing that the DGRB due to the CR interac- 
tions in clusters should be a small fraction of the total diffuse flux observed 
by EGRET. This conclusion was reached by the previous authors under the 
hypothesis that a large fraction of the baryons in the universe is contained 
inside clusters of galaxies (BBP considered that clusters are a fair sample of 
the baryons in the universe [17], [18], [19]) assumed to have a homogeneous 
inner distribution of gas, n e = const (here n e is the IGM electron number 
density). Because of these assumptions, their results depend only on overall 
cosmological parameters like the baryon fraction in the universe Qb, and on 
the cluster size. 

Dar & Shaviv [12] also predicted the gamma ray fluxes from a few nearby 
clusters (Coma, Perseus and Virgo): for these three clusters they found 7-ray 
fluxes in the range F 7 (> 100 MeV) ps 5 4- 20 ■ 10 _8 cto _2 s _1 . In particular, 
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the value which they predicted for A1656 (Coma), F 1 (> 100 MeV) ~ 5 • 
10~ 8 cm _2 s _1 , is close to - or slightly higher than - the upper limits given by 
EGRET for this source. Similar results were obtained for these clusters by 
Ensslin et al. [20] assuming a population of CRs from radio sources located 
within galaxy clusters in almost equipartition with the IGM thermal energy. 

We stress here that in all the previous works a uniform IGM density profile 
was assumed. Moreover, the cluster population was not assumed to evolve with 
cosmic time, and the same working hypothesis of no-evolution was assumed 
for the IGM content of each cluster. 

However, X-ray studies of galaxy clusters, have shown that these cosmic 
structures are indeed well structured, having a gas density profile n(r) oc 
[1 + (r/r c ) 2 ]~ 3/3 / 2 , with core radii r c « 0.1 4 0.3 hr l Mpc and (3 « 0.6 4 0.8 
(see e.g. [21]; see also [15] and references therein). Beside this, the IGM is in- 
deed evolving as indicated by its sensitive metal enrichment, Fe/H ~ 0.24-0.5 
(in solar units, see e.g. [22], [23]), shown even for the brighter clusters observ- 
able at redshifts z ~ 0.5 [24]. Nonetheless, there is also an increasing debate 
on the possible evolution of the X-ray luminosity function observed out to 
z < 0.5 with the EINSTEIN [25] [26] and ROSAT satellites [27] [28] and on 
the possible evolution of the cluster temperature function [29] [30] [31]. If an 
evolution is present in the cluster population this can be, in fact, understood 
as a result of two competing effects: 

i) a luminosity evolution, where the cluster X-ray luminosity, L oc n 

2 T 1/2 R3 

(mainly due to thermal bremsstrahlung) , changes with redshift due to varia- 
tions in the gas mass density, n oc f g p c i (where the cluster gas mass is taken 
to be a fraction f g = M gas /M of the total cluster mass), and/or changes in 
the IGM temperature T at fixed mass, M oc p c iR 3 (here p d is the cluster total 
mass density); 

ii) a change in redshift of the number density, N(M,z), (usually referred to 
as mass function, hereafter MF) of clusters that are found to be collapsed (or 
virialized) in the mass range M, M + dM at redshift z. 

Detailed studies of cluster evolution in X-rays (see e.g. [32], [33], [29]) con- 
sidered in fact that a combination of the previous mechanisms is responsible 
for the actual cluster evolution when they fit the available data (see [32] for a 
detailed discussion). 

In this paper we predict the amount of high energy, non-thermal, gamma-ray 
emission from galaxy clusters using detailed modelling of the realistic cluster 
structure, as well as viable modelling for the evolution of the IGM and of the 
cluster MF. Based on these phenomenological cluster models, we predict the 
amount of DGRB that can be produced in the viable cosmological models: here 
we consider flat and low-density (open or vacuum-dominated) CDM models 
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as well as mixed Dark Matter models with a fraction Q v pa 0.3 of the total 
density of the universe in form of massive neutrinos. We use h = ifo/100 km 
s _1 Mpc -1 throughout the paper unless otherwise specified. 

The plan of the paper is the following. In Sect. 2 we briefly summarize the 
cluster formation hystory in hierarchical scenarios for structure formation. In 
Sect. 3 we describe a model for the production of diffuse gamma-ray emission 
due to the interaction of CR's with the target protons present in the extended, 
diffuse IGM. We consider in Sect. 4 different CR sources that can be found in 
connection with galaxy clusters. We discuss in Sect. 5 the correlation between 
the extended gamma-ray emission and the much better known thermal X-ray 
emission coming from the IGM. Based on these properties, we construct a 
list of predicted 7-ray fluxes for a compilation of X-ray clusters with detailed 
informations on their IGM structure, IGM temperatures and X-ray fluxes. In 
Sect. 6 we present predictions for the amount of DGRB produced by galaxy 
clusters in different cosmological scenarios. We briefly discuss in Sect. 7 the 
extended neutrino fluxes emerging from these objects and their contribution 
to a possibly detectable diffuse neutrino background (hereafter DNB). Finally, 
in Sect. 8 we discuss our results in the light of the current limits obtained from 
EGRET and in the light of the future experiments for gamma-ray and neutrino 
astronomy. 



2 Cluster formation: theory 

Clusters of galaxies originate from small density perturbations at early times 
[34] that grow under the action of their own gravitational instability. In the 
viable cosmological scenarios, cluster evolution can be followed adopting a 
simple spherical collapse picture according to which a homogeneous, spherical 
perturbation detaches from the Hubble flow at time t m given by the equation, 
t m = [37r/32Gpfe(t m )] 1 / 2 , collapses at time t c ~ 2t m , and virializes at time 
t v ~ 3t m [here pb{t m ) is the cosmological background density at the time 
t m of the maximum expansion of the perturbation, see [69] for details]. The 
relative density contrast of such a perturbation at the initial redshift Zj, 8 i>v , 
depends on the cosmological model. Under the assumption of linear growth, 
the density contrast at t v is 8 V = 5i >v V(t v )/V(ti), where T>(t v ) is the linear 
growth factor in the chosen cosmology (see [35] for details). For flo 1, 5 V 
tends to the standard value 5 V ~ 2.2, independent of t v . The actual, non- 
linear density contrast, A = p/p b (where p is the perturbation density and 
Pb is the background density at the time of virialization) of a cluster that 
virializes at redshift z in a Qq < 1 cosmological model writes: A(f2o, z) = 
187r 2 /[r2o(^o^) 2 (l + <2) 3 ] [ m flat, vacuum dominated low-density models A has 
not an analytical expression (see, e.g. [35] and references therein)]. It tends to 
the standard value ps 400, found in a flat (Q = 1) universe. 
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Detailed imaging of galaxy clusters in the X-rays (see [15]) revealed that the 
IGM is concentrated in the cluster central region (the core) and its density de- 
creases rapidly at large distances from the core. Thus, following Colafrancesco 
et al. [35], we relax the assumption of uniformity, n = const, by considering 
a 3-D gas density profile: 



n(r) = n c 



1 + 



-3/3/2 



(3) 



where n c is the central electron number density and r c is a core radius. 

There are various open issues pertaining to the formation of hot gaseous cores 
of galaxy clusters. For our purposes here we can adopt the following simpli- 
fied approach (see also [35]): shortly after a cluster forms and virializes, a 
gaseous core forms (probably as a result of tidal galactic interactions and 
other gas stripping processes) with the hot gas in local hydrostatic equi- 
librium in the potential wells of the cluster. Analytical studies of the clus- 
ter self-similar collapse [36] and numerical simulations [37] [44] indicate that 
the gas density profile scales following the total matter density of the clus- 
ter. Thus, the central gas density, n c , is related to the total density at the 
cluster centre, p c , through n c = f g (p c /m p )2/(l + X), where m p is the pro- 
ton mass and X = 0.69 is the cosmic Hydrogen mass fraction. The mass 
within the outer radius, taken as R = qr c , is M(q,f3) = 3M c u(q, (3), where 
^{q, P) = Jo dxx 2 (l + x 2 )~ 3f3 ^ 2 . Here M c = [An /3)r 3 p c . Because of the assumed 
profile, the ratio between the central and mean mass density of the cluster is 
p c /(p b A) = q 3 /3u(q, (3). Assuming that the cluster collapse is self-similar, we 
can infer the mass and redshift dependence of the core radius r c from the 
equivalence: r c = R v /q = i(3M/47rp 6 A) 1 / 3 (l + This gives: 

„ , l^/i^Mpc 
r c (n , M, z) = — 



M A(1,U) ^° 1 

io 15 !!- 1 ' KKlnZz). i + z 



■ (4) 



Hereafter we fix q = 10 to recover r c values consistent with the observations 
[47] [21]. 

Under the standard assumption of the IGM in hydrostatic equilibrium with 
the potential well of a spherically-symmetric virialized cluster, the relation 
between IGM temperature and cluster mass is easily obtained by applying the 
virial theorem: T = —pm p U/(3kM), where p = 0.62 is the mean molecular 
weight (corresponding to a hydrogen mass fraction of 0.69), k is the Boltzmann 
constant and U is the cluster potential energy. If the cluster is assumed to be 
uniform, U = -(3/5)GM 2 /R v and T = = (l/5)(pm p /k)GM/ R v , where 
R v = [3M/(47rpfeA)] 1 / 3 /(l + z) is the cluster virial radius. In the case of the 
density profile in eq. (3), U = -(GM 2 /r c )^/u 2 and T = 5#T^/ '(3lu 2 ) , where 
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rj} = f«ydy/(l+y 2 ) 3m Jo t 2 dt/(l + t 2 ) 3m - Thus the cluster IGM temperature 
is 



T = 5.8 keV (1 + z) 



M \ 2 / 3 



HoAMi 1/3 

A(1,0) J ' {b) 



where the normalization constant, T = 5.8 keV, is obtained for q = 10 and 
(5 = 2/3. For (5 — 1 we get T ~ 8.2 keV, where T is the temperature of a 
1O 15 M /i _1 cluster at z = in a f2 = 1 universe. 

Although it is expected that the IGM mass fraction, f g , depends on z and M, 
little is currently known on the exact form of these dependences. We adopt the 
simple parametrization (described in detail in Colafrancesco & Vittorio [32]) 
which is suggested by models of the IGM evolution driven by entropy variation 
in the cluster cores [38] [39] [40] and/or shock compression and heating [41] 
[42] [43]: 

The normalization to f 9t0 ~ 0.1, is based on a local, rich cluster sample. Values 
of rj ~ 0.2 -f- 0.9 and s ~ 0.5 -r- 2 are consistent with the available data. 

The previous information about n(r) [see eq.(3)], T [see eq.(5)] and the clus- 
ter extension, R = qr c , allow to predict both the gamma-ray and the X-ray 
emissivities which will be discussed in the following sections. 



3 Diffuse High Energy Emission from Clusters of Galaxies 



As recently discussed by BBP [13] and Volk et al. [14], most of the cosmic rays 
produced in clusters of galaxies remain confined within the cluster potential 
wells and produce high energy gamma rays and neutrinos by interactions with 
the intracluster baryonic gas (the IGM). The confinement of these CR's in 
the intracluster space is a crucial mechanism for maximizing the efficiency of 
the 7 (or neutrino) emission and is strictly related to the value and to the 
configuration of the magnetic field in clusters. This determines the diffusion 
coefficient, given by: 

D(p} = r L(p) ifk5^m- (7) 

Here P(k) is the power spectrum of the fluctuations in the magnetic field and 
tl(p) = pc/(eB) is the Larmor radius of a particle with electric charge e and 
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momentum p. In this paper we assume that the fluctuations of the magnetic 
field have a Kolmogorov spectrum: 



/ k \ ~ 5/3 

P(k) = P(ko) ( -j- J , (8) 



where ko = 1 /do, and do is the smallest spatial scale at which the magnetic 
field achieves homogeneity. The spectrum in eq.(8) is normalized as: 

+00 

J P(k)dk « B 2 , (9) 

ko 



which implies, 



- 1 ~>2 _ 2 J r>2 



P(fco) = --^ = -rf 5 2 . (10) 



Finally, the diffusion coefficient for relativistic particles is obtained from eqs. 
(7) and (10): 

D CR (E) = \cdT{eB)-^E^ , (11) 



where E = pc is the energy of the relativistic particles. 

Measurements of the overall magnetic field in single galaxy clusters yield values 
in the range ~ 1 4- 10 /iG ([49] [50] [51] [52] [53]), while a statistical analysis 
of a sample of clusters (see e.g. [54]) yielded typical, average values B ~ 1/j.G 
on homogeneity scales do ^ 20 kpc. This size is compatible with the model 
of Jaffe [55] for the origin of the magnetic field in clusters: according to this 
model, the turbulence of the IGM produced by the large scale motions of the 
galaxies is responsible for the value and the homogeneity scale of the magnetic 
field. 

Thus, from eq. (11) we obtain: 

D CR {E) w 2.3 x 10 26 E{eV) 1/3 B^ /3 cm 2 / s . (12) 



The diffusion time of particles with energy E at distance r from the CR source 
reads: 



r 2 



T = 



QDUEj * 6 ' 9 X ™ 21 r 2 M PC E(eV)-^s, (13) 
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where the factor 6 in the previous expression for r comes from the solution 
of the diffusion equation in spherical geometry (here we assumed B^q = 1). 
This implies a confinement of cosmic rays, on typical scale of R c i ~ 2 Mpc, if 

E £ 4.2 x 10 14 /T 3 eV . (14) 



These cosmic rays have diffusion times larger than the age of universe, for 
which we adopted the value to = 2.06 x lO 17 ^" 1 s. 

More difficult is to extrapolate, from the average confinement ability of the 
clusters, the confinement efficiency of the cluster core: the difficulty comes 
from the small size of this region (r c ~ 0.2^-0.5 Mpc) where a detailed knowl- 
edge of the structure of the magnetic field is needed (the typical separation 
between galaxies in clusters cores can be even smaller that the size of the 
largest galaxies often observed in the cores). 

The confinement of cosmic rays inside clusters and the relevance of this process 
for the production of high energy radiation was already pointed out by BBP. 
In that work, however, the true density profile of clusters was not considered, 
and only an average IGM density entered the calculation. Indeed, for realistic 
IGM density profiles, like that in eq. (3), a large fraction of cosmic rays could 
be confined inside the dense cluster cores, thus determining a sensitive increase 
in the interaction rates. 

Following the approach of BBP, and assuming a power-law spectrum of CR 
produced by a source located at the center of the cluster, the number of i- 
secondaries (i= 7, v) produced per unit time and unit volume, at energy E 
and at distance r from the CR source in a cluster, is: 

qi(E,r) = Yi(r pp n 2 (r)c 



n p {ti,r) 



n(r) 



(15) 



where a pp ~ 



3.2 x 10 cm , Yi are the yields for the i-secondary production 
n(r) is the IGM density profile given by eq.(3) and the produced CR 
proton density, n p (E,r), is determined from the diffusion equation as: 



n p (E,r) 



" Qp(E) ■ 
-Dcr(E). 



(16) 



In this last equation, Q P (E) is the emitted spectrum of CR from a source as- 
sumed to be located in the cluster core. Thus, the total number of i-secondaries 
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produced per unit time at energy E reads: 



R 



Qi{E) = 4tt J drr 2 qi(E,r) = YiCr pp n c 



Q P (E) 



D C r(E) 



r 2 c aq,(3,E,z;M) , (17) 



where 

Rdiff/rc 

((q, P, E, z; M) — J dx x(l + x 2 )~ z<il2 . (18) 
o 



We note that the maximum length over which CR diffuse is Rdiff = (-d — 
\j6Dcii(E)t(z, Q ), at each epoch t(z,Q ). This yields 



2-10 1 V/ V10 29 cm 



For a uniform cluster we reproduce the results of BBP. In the present approach, 
instead, the inclusion of a realistic density profile for the IGM produces a 
change from a pure power-law CR spectra, Qi(E) oc E~ l9 [if Q P (E) oc E~ la 
is assumed], due to the presence of an energy dependent term in the function 
C oc [1 + (^d/?"c) 2 ]~ 3/ ^ 2+1 in eq.(18). However, at the energies for which the CR 
confinement is mostly effective, <^ r 2 , expanding eq.(17) in power series, 
yields: 

Qi(E) w ZY l( r pp n c ct{z)Q p {E) . (20) 



In this limit, the total number of i-secondaries produced per unit time is 
independent on the cluster size and depends only on the IGM density, n c at 
the cluster centre. 

Note that the term Q P (E) can take into account the possible evolution of the 
luminosity of the CR source as well as the possible evolution of the number 
of CR sources that are present in the central regions of the clusters. Another 
important source of evolution in the 7 (or neutrino) diffuse emission from 
galaxy clusters comes from the possible evolution of the IGM content (see 
Sect. 2). We will discuss these points in the following. 
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4 Sources of Cosmic Rays in Clusters of Galaxies 

In the previous sections we discussed the confinement of cosmic rays in galaxy 
clusters, irrespective of the way they are produced. In this section we shall 
deserve more attention to the possible CR sources and we will estimate their 
contribution in terms of CR luminosity. 

4-1 Normal galaxies 

The most natural candidates as sources of CR in clusters are normal galaxies. 
If we assume that our own galaxy is a typical one, the CR luminosity can be 
estimated to be Lcr — 3 x 10 40 erg/ s [56]. The central number of galaxies Q 
iVo (which is a measure of the cluster richness) is found to be correlated with 
the cluster temperature and it scales as N oc T as [22]. Thus the gamma ray 
luminosity associated with normal galaxies in clusters scales as Lcr oc N L p oc 
T 8 L p , and for a Coma-like cluster with T = 8.3 keV, we find Lcr ~ 2 • 10 42 
erg/s. In any case, even for the richest and hottest clusters observed (A2163 
with T = 13.9 keV) the CR luminosity associated to normal galaxies is found 
to be L C r ^ 10 43 erg/s. 

4-2 Active galaxies 

As discussed by BBP, on statistical grounds we can expect ~ 1 active galaxy 
per cluster, where the term active galaxy is used here to indicate several 
classes of luminous objects, such as AGN, radiogalaxies and cD galaxies. 
The mechanisms of CR production in active galaxies have been studied by 
many authors (see e.g. [57] and references therein) and typical luminosities of 
Lcr ~ 10 44 erg/s can be achieved. 

In a recent paper, Ensslin et al. [20] proposed that jets of powerful radiogalax- 
ies in clusters can establish a sort of equipartition between the thermal energy 
of the IGM and the cosmic ray energy density in clusters, due to particle 
acceleration in the jets. The previous authors calculated the gamma ray flux 
from single clusters due to the CR interactions in the case mentioned above 
and predicted a gamma ray flux of F^(> lOOMeV) ~ 6 x 10~ 8 cmr 2 s~ x for 
A1656 (Coma), a little bit in excess of the EGRET experimental limit. Sim- 
ilar predictions are given for other clusters (A426, Ophiucus, M87- Virgo) for 
which they found F 7 (> WOMeV) in the range m 3 12 X 10~ 8 cm~ 2 s _1 . 

1 The central number density of galaxies Nq is related to the total number of 
galaxies in a cluster by the expression N g = 4TrNor^uj(q, (3). 
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Given all these estimates, a CR luminosity Lqr ~ 10 44 erg/s can be expected 
from active galaxy interaction with the IGM of local clusters. At higher red- 
shifts, the AGN luminosity (or the CR density) is expected to increase roughly 
as ~ (1 + z) 2 while the IGM density (the target proton density) is expected 
to decrease roughly as ~ (1 + z)^ 1 - 2 . The two evolutionary effects tend to 
balance out and the net result is a moderate increase, if any, of the effect due 
to distant, AGN populated clusters. 

4-3 Cluster accretion shocks 

Another important acceleration mechanism for CR's is provided by accretion 
shocks formed in and around clusters of galaxies as a result of their collapse and 
virialization [36] [58] [59] [60]. These shocks have been recognized as possible 
sources of CR acceleration through the first order Fermi mechanism (see e.g. 



We shall estimate here the CR luminosity by adopting two models for the 
accretion of matter around galaxy clusters. In the first model we refer to the 
analytical calculations of Bertschinger [36], in which a self-similar analytical 
solution of the hydrodynamic equations is found. The solution, concerning 
the velocity and density of the accreting gas, is written in terms of the di- 
mensionless variable, A = r/R ta , where r is the distance from the center of 
the cluster (assumed to have spherical symmetry) and R ta is the turn- around 
radius, given by 



(here M is the cluster gravitational mass). In these calculations [36] the sta- 
tionary shock is located at \ s h ~ 0.347. At the shock location, the gas ve- 
locity is found to be v(\ s h) ~ 1.5R ta /t , while the density of the gas is 
Pgas(^sh) ~ 2.2 x 10~ 29 g/cm 3 . This last number has been obtained assuming 
Pgasij') = 0.1p(r), where p is the total mass density of the cluster (dark mat- 
ter and baryonic gas have been assumed to have the same radial distribution). 
Thus, the total energy per unit time available at the shock can be written as: 



(we use here h = 0.75). The typical efficiency of CR acceleration at the shock 
is of order ~ 0.1 [56], so that Lqr ~ (2-^3) x 10 44 erg/s is obtained, which is 
comparable with the CR luminosities derived in the previous Section 4.2 for 
the case of active galaxies within galaxy clusters. 



[13] [59]). 




(21) 



L = \p g as{Rsh)v{R sh )HiiR 2 sh v{R sh ) » 2.4 x 10 



erg/s 



(_Rta_\ 5 (22) 

\5MpcJ K J 



12 



An alternative way to estimate Lqr for the accretion shock case is to write it 
through the matter accretion rate, M, at the shock. Numerical simulations of 
the large scale structure formation ([61] [62]) give larger values for the shock 
size, of order ~ 5 h^Mpc. If we again assume that a fraction ~ 0.1 of the 
gravitational energy at the shock is converted in cosmic rays accelerated at 
the shock, we can write: 

GMM , iA 

L C r ~ 0.1— « (1.3 - 2.5) x 10 44 erg/s. (23) 

R s h 



In this numerical estimate, we assumed that the cluster accretion rate at the 
shock is given by the time averaged value 

M^<M>= ~ 5 x 10 29 g/s, (24) 



where M gas is the mass of the cluster in the form of IGM, assumed here to 
be ~ 10% of the total mass. Note that this last estimate is probably a lower 
value for the actual IGM mass in rich clusters (see e.g. [18] [63]). 

To summarize, the CR luminosity obtained in the two accretion shock models 
yield comparable results for the local cluster CR luminosity, Lqr ~ a few 
10 44 erg/s. At higher redshifts, the CR luminosity from cluster shocks should 
tend to decrease, if f g does. 

On the basis of the previous considerations we are justified in assuming an 
average CR luminosity of Lcr = 10 44 erg/s for a typical, rich cluster with 
M = Mi 5 (here M± 5 = 10 15 h _1 M ) at z — 0, either due to active galaxies 
and/or to accretion shocks around clusters. 



5 The correlation between X-ray and Gamma-ray emission 



From eq.(17) we can now calculate gamma-ray fluxes and luminosities for 
clusters with mass M at redshift z. In Fig.l we plot the cluster 7-ray luminosity 
against the cluster temperature T expected in different, viable cosmo logical 
scenarios (see Fig.l caption for details). In our specific predictions we find 
that the cluster 7-ray luminosity L 7 [or equivalently Qi as written in eq.(17)] 
scales with the cluster properties (at fixed redshift z) as: 



Ly OC 



QpjE) 
D C r(E) 



n r 2 c ( , 
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that yields the scaling: 



L 1 oc 



• Q P (E) 
-Dcr(E) 



(25) 



We assumed here a virial equilibrium for the IGM which yields, T oc GM/R. 
In our specific model for the IGM we also consider a distribution of values for 
f g = M gas /M oc M v (with rj ~ 0.2 4- 0.9) from groups to galaxy clusters in 
order to recover phenomenologically the general trend indicated by the most 
recent data [63] [64] [48] . The function ( decreases with increasing cluster mass 
(temperature) and hence the actual scaling of L 7 with the cluster temperature 
reads L 1 ~ T 5 that is weaker with respect to the scaling of the previous 
equation (see also Fig.l). The shape of the L 7 (T) curves also depends slightly 
on the considered cosmological model, with low Q values providing higher 
L y at fixed T (see Fig.l). This mild cosmological dependence is caused by 
the sensitivity of the cluster parameters (mainly the core radius r c ) to the 
underlying cosmology: in fact, r c increases with decreasing Q [see eq.(4)]. 

In Fig. 2 we plot the expected 7-ray fluxes at E > 100 MeV for clusters with 
X-ray flux Fx{2 — lOkeV) in the energy band 2 — 10 keV. We show this plot for 
clusters at z — 0.023, the redshift of Coma. For this cluster we have an upper 
limit on F 7 (> 100 MeV) measured by EGRET (as indicated in Fig.2). We 
plot the fluxes expected in the cosmological models considered in Fig.l. The 
7-ray luminosities shown in Fig.l correspond to fluxes F 7 (> 100 MeV) & 10~ 9 
photons s" 1 cm" 2 for clusters at z ~ 0.1, and F 7 (> 100 MeV) in the range 
7 • 10~ 9 -i- 1.5 ■ 10~ 8 photons s" 1 cm -2 for Coma-like clusters (at z ~ 0.023) 
in different cosmological models, as shown in Fig.2. These fluxes are a factor 
~ 2 -f 5 smaller than the EGRET upper limit for Coma; such a limit was 
saturated by the predictions of Dar and Shaviv [12] and Ensslin et al. [20]. 
Smaller fluxes are obtained for flat (Qq = 1) cosmologies. The differences in 
L 7 expected in different cosmological models (see Fig.l) are amplified when 
one compares the relative fluxes, F 7 = L^/4nd 2 L , by the fl dependence of the 
luminosity distance di- At redshifts z ^ 0.3 the predicted -F 7 (> 100 MeV) 
for clusters fall below values ~ 10~ 10 photons s _1 cm" 2 , and become hardly 
observable from any present and/or planned gamma-ray mission for the next 
coming years (INTEGRAL, GLAST, AMS). 

From the previous results, we found a correlation between the X-ray luminos- 
ity of galaxy clusters and their gamma (or neutrino) luminosity due to CR 
interaction in the intracluster gas, as shown in Fig.2. Such a correlation is 
indeed expected because the protons in the IGM are both the target for the 
CR interactions and the responsible for the thermal bremsstrahlung radiation 
which provides the bulk of the cluster X-ray luminosity. In fact, the total X-ray 
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luminosity from a cluster with an IGM temperature T can be written as: 

R 

L x = A(T) (kT) x ' 2 J dr4nr 2 n 2 (r) . (26) 
o 

Using eq. (3) for the IGM density profile, eq.(26) writes as: 

L x = Aixny c A{T){kTfl 2 l{q^) , (27) 

where 

Q 

l(q,/3) = ( dxx 2 {l + x 2 )-^ . (28) 



The function A(T) in eqs. (26-27) contains the contribution of the Gaunt 
factors to the frequency integrated X-ray spectrum, and writes as: 

2 5 7tp 6 / 27r \ 1 / 2 r 



where y = hu/kT, e and m e are the electon charge and mass, respectively, 
and gff is the Gaunt factor. 

Using eqs. (8) and (16) we can derive the expression for the ratio of the 
the gamma (or neutrino) flux coming from a single galaxy cluster and its 
corresponding X-ray flux, Fx- 

Fy tU Y^ v <7ppC / 1 \ 1 I Qp(E^) 



(n c r c ) (kTV/A Dc R (E) C \ ■ (30) 



F x 4nA(T)l(q,[3) \n c rj (kT) 1 / 2 [D C r(E) 

Thus, from the previous eq.(30) it is possible to estimate the cluster 7 (or 
neutrino) flux given a measurement of the cluster X-ray flux Fx, the central 
IGM density n C) the cluster core radius r c and the cluster IGM temperature 
T: 

F„ = F x .C(E,/3,z) 1 -= (31) 



Here the quantity C(E,f3,z) reads: 

\ 

4^A(T)I(q,P) \ D CR (E) 



C(E,P,z) = - Y A Z^f m { ^ E L c) ■ (32) 
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We will use here the previous results and the available data for Fx,n c ,r C} (3 
and T [21] [65] [67] to derive a list of estimated cluster fluxes in the high 
energy 7 and neutrino energy regions (see Table 1). The clusters in the sample 
listed in Tab.l are those which have well measured values of (3,n c ,r c ,T and 
F x and are selected from the homogeneous sample of Jones and Forman [66] 
which is limited to redshifts z ^ 0.1. In Fig. 3 we plot these data in the -F 7 (> 
lOOMeV) — T plane, together with the predicted curves, F 7 (> 100MeV)(T), 
expected in a CDM+A model with fl = 0.4 and h = 0.6 at redshifts z = 0.023 
(solid line) and at z — 0.072 (dashed line). The theoretical predictions at 
the minimum and maximun redshifts of the clusters in the Jones & Forman 
sample encompass the data points. The 7-ray fluxes predicted from eq. (30-31) 
for our sample of clusters at z 0.1 span over the range F 1 (> lOOMeV) ~ 
5 • 10~ 10 -=-8.5 • 10~ 9 photons s _1 cm -2 and are distributed (with some intrinsic 
dispersion) along the theoretical F 7 (> lOOMeV) — T curves predicted at each 
redshift (see Fig. 3) according to eq.(31). Note that the region at T 2 keV is 
populated mainly by galaxy groups for which the determinations of the IGM 
temperatures are more difficult and uncertain due to their small amount of 
diffuse IGM (see [45] [46] for a discussion). In the following, we will consider 
mainly galaxy clusters with T ^ 2 keV. 



In Fig. 4 we show how the clusters listed in Tab.l distribute according to 
their predicted 7-ray flux. The distribution of the cluster 7-ray fluxes peaks 
at F 7 (> lOOMeV) ~ 8 • 10~ 10 3 • 10~ 9 photons s^ 1 cm" 2 . Our sample of 
clusters was extracted from an X-ray flux limited sample [66] and thus it can 
be considered as representative of the cluster population even in the 7-rays 
because of the existing correlation Fj — F x of eq.(32). From Fig. 4 we note that 
the optimal observational strategy to detect large samples of galaxy clusters 
in the 7-rays is to achieve sensitivities ^ 5 • 10~ 10 photons s _1 cm~ 2 : such a 
sensitivity level should be at hand with the next generation 7-ray telescopes 
(INTEGRAL, GLAST). 



6 The contribution of galaxy clusters to the DGRB 



The diffuse flux of 7 (or neutrino) emission at energy Ei that is received from 
a cluster at redshift z is given by 

F.AEr) = ^ , (33) 
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where the luminosity distance 



z 

d L (z,Q Q ) = ^-(l + z) fdz' 
H J 



-1-1/2 

(1 + z'ftt + l-n , (34) 



allows for the flux-luminosity conversion: 

L = AnFd 2 L (z, Q ). (35) 

The total flux of 7-rays (or neutrino) due to the superposition of a population 
of evolving clusters is then: 

r / p a Z J x dV{z,Q ) 7 Qi\Bh{l + zlz\ 
I„{E)= j ~^—dz j dMN(M ) z) ^ , (36) 



where dV(z,flo)/dz is the cosmological volume element per unit steradian. 

The space density of clusters at any mass M and redshift z is given by the mass 
function, N(M,z), (hereafter MF) usually derived by the Press & Schechter 
[68] theory as: 



We remind the reader that pb is the co moving background density of the 
universe, M is the total cluster mass and 5 V is the linear density contrast of 
a perturbation that virializes at redshift z (see Sect. 2). The variance a of the 
(linear) density fluctuation field at the scale R = (3M/A-Kp b ) l l z and redshift z 
is given by the standard relation: 

a , {RiZ ^m^j edkm ^y (38) 



[69], where V is the growth factor of linear density fluctuations in a given 
cosmology and ji(x) is a spherical Bessel function. If we normalize the fluctu- 
ation spectrum, P(k), by requiring a (8 h- 1 Mpc,z = 0) = ft -1 , then the MF 
depends uniquely on the product b5 v , for a given cosmological model. 

Assuming a power-law power spectrum for the initial fluctuation field, V(k) = 
Ak n , the MF attains the following analytical expression: 



N(M, z) 




Pb n + 3 S v b ( M 



6 Viz) \Mn 



cxp 



1 by rM\ 

2V\z) \MJ 



2n 



, (39) 
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where a = (n+3)/6 and Mq is the mass contained in a 8 h~ x Mpc radius sphere 
at z — 0. The MF we will use here is normalized to the observed abundance 
of clusters at z ~ by fitting the parameter S v b (see [32] [29]). This fitting 
procedure is done by using the available data on the temperature function [67] 
and/or those on the X-ray luminosity function [70] [27]. The different data sets 
yield approximately the same best fit values for 8 v b in the cosmological models 
we considered. Specifically, CMV [29] found 8 v b = 2.5 for a standard CDM 
model and 5 v b = 1.6 for a low density, vacuum dominated CDM model with 
Vl = 0.4 and f2 A = 0.6, which is the best fitting CDM model. For this class 
of CDM models, CMV found 5 v b « 2.5fi°- 49 , or equivalent^ a s « 0.67fio° 49 
(if one assumes S v = 1.68). 

In the case of power-law power spectra and for Q = 1, eqs. (34), (36) and (39) 
give the following analytical result for the cluster contribution to the DGRB: 

x J d z(l + z)-^-^(J dyy-^exp \~[^J V 

Vmin V 



where y = M/M , x — r c(l + ^)/M 1 / 3 has been obtained from eq. (4), ( 
is defined in eq. (18) and e is the power index of the energy dependence of 
D C r(E): Dcr(E) oc E € , where we found e = 1/3 (see eq.ll). 

A numerical calculation of the integrals in eq.(40) yields the immediate result 
that the main contribution to the 7 (or neutrino) diffuse flux comes from low- 
z clusters, and actually from sources at z ^ 0.3. This result is confirmed by 
our numerical calculations of the redshift distribution, N(z), of 7-ray emitting 
clusters found in the different cosmological models we considered in this paper. 
In Fig. 5 we show the redshift distributions, N(z), as a function of z, of 7-ray 
clusters expected in a flat (Q Q — 1, h — 0.5) CDM model in the case of an IGM 
evolving as f g oc M°- 2 (l + z)~ 1 ' 8 (solid histogram) and in the case of absence of 
IGM evolution with redshift (dotted histogram). We considered clusters with 
fluxes .F 7 (> lOOMeV) > 10~ 10 photons s -1 cm~ 2 with 7 9 = 2.1. The presence of 
an IGM evolution modifies the cluster distribution mainly at z <^ 0.2. Models 
with no IGM evolution (i.e. with s = 0) predict ~ 10 times more clusters 
with F 7 (> lOOMeV) Z 10~ 1C ^photons s" 1 cm' 2 at z k, 0.4 with respect to the 
case of a strong evolving IGM model with s = 1.8. However, such a difference 
tends to reduce to ~ 30% at z ^ 0.2 and becomes negligible at z ^0.1. As 
the bulk ( k 90%) of the bright 7-ray clusters is located at z ^ 0.15, the IGM 
evolution should play a little role on the cluster contribution to the DGRB. 
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The value of z max in eq.(40) is denned by two criteria: i) the largest cutoff 
in the redshift distribution of the cluster population; ii) the maximum energy 
of the i-secondaries. In fact, for E ^ E max , the energy of production of the 
i-secondaries is E(l + z) ^ E max (that is the maximum energy of the original 
CRs). If E > E maX) we must realize that the spectrum of the i-secondaries is 
changed. 

Note that in these reasonings, the diffusion coefficient D C r(E) depends on 
the typical homogeneity scale d of the magnetic field in clusters (see Section 
3). There are reasonable arguments to expect that this scale could change 
with redshift. However, we assume here that the only dependence with z of 
the diffusion coefficient Dcr is due to its dependence on the cosmic time (see 
eqs.11-13). 

In Fig. 6 we show the cluster contribution to the DGRB, 7 7 , as a function of the 
7-ray energy E, in the various cosmological models we consider in this paper 
(see Fig. 6 caption for details). The expected contribution of galaxy clusters 
to the high energy DGRB [evaluated according to eqs. (36-40)] is found to be 
^ 5% of the EGRET value [1] for all the considered cosmological models (see 
Fig. 6). This result depends on different issues: 

i) we choose Lqr = 10 44 erg/s for a M15 cluster in each cosmological model. 
The value of J 7 depends linearly on the choice of Lcr so that any variation of 
this parameter reflects in a analogous variation of the predicted contribution 
to the DGRB; 

ii) we normalize the local abundance of clusters to that observed in the X- 
rays. This estimate is the most precise and robust at the moment. Such a 
normalization yields a values of 8 v b for each cosmological model, so that the 
evolution of galaxy clusters is completely determined; 

Hi) the distribution of the cluster IGM is self-consistently derived from the 
gravitational instability picture for the formation of galaxy clusters. The level 
of 7 7 depends linearly on the value of f g o(Mis, z — 0) for which we choose the 
value 0.1. Changing this parameter in the viable range, f g0 ~ 0.05 — 0.3 yield 
similar variations in J 7 ; 

iv) only \ow-z clusters effectively contribute to the integral in eq.(36) as shown 
by the redshift distribution of clusters with F 7 > 10 -10 cm^s' 1 (see Fig. 5). 
Nonetheless, there is a mild dependence of the results from the considered 
cosmological model. The MDM models show the smallest value of the DGRB 
contributed by clusters due to the very recent formation epochs predicted for 
clusters in this model (see Fig. 6). CDM models (flat and low density) show 
approximately the same level of DGRB: this is due to the normalization of 
their spectral amplitude to the observed abundance of nearby clusters observed 
in the X-rays [29]. 
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In Fig. 7 we show the quantity J 7 as a function of the 7-ray energy E, expected 
in a CDM+A model (fl = 0.4, h = 0.6) for the cases of the maximum rate 
of IGM evolution allowed by the present data, f g oc M 02 (l + z)~ 2 - 2 (solid 
line) and for the case f g = const (dashed line). We note that the amount of 
DGRB contributed by galaxy clusters in different cosmological models depends 
slightly of the IGM evolution (see Fig. 7) as the clusters effectively contributing 
to the DGRB are those with z ^ 0.2. At these redshifts, the maximum change 
of f g (at fixed mass) is of the order of ~ 20% and this corresponds to a 
variation of I 1 of the same order (see eqs.18, 36 and 40). 

Variations in the parameters of the cluster structure yield variations ^ 5% on 
the final result. Changing /3 by 20% yields variations in ( of ~ 8% that give 
similar changes on 7 7 . Variations in the primordial index n of the perturbation 
spectra yield variations of order ~ (n + 3)/6 on 7 7 . 

Thus, while the level of the DGRB predicted for a specific choice of the cluster 
evolution models are at the level of 0.5-^0.8% of the EGRET diffuse flux, more 
conservative estimates of J 7 require to consider the theoretical uncertainties 
in the modelling of cluster structure and evolution. 

Theoretical uncertainties in the description of the relevant quantities of cluster 
structure yield variations of the results contained within ^ a few %. Let us 
be more specific. 

Our results are based on the theoretical description of cluster evolution through 
the spherical top-hat model and the PS mass function. Despite the inherent 
uncertainties in the detailed models for cluster collapse, the final results de- 
pend on the combination bS v that is fixed by the fitting procedure to the local 
cluster abundance with uncertainties ~ 5 -j- 8% [29]. Moreover, our results on 
7 7jl/ are weakly sensitive to the lower cutoff, M min , in the MF (see eqs 35-39). 
This reduces the effects of possible changes of the slope of the MF in the efforts 
to go beyond the simple, first-order theory for the origin of the universal MF 
(see [48] for a discussion). 

Uncertainties in the description of the IGM evolution (both its normaliza- 
tion, f g0 ~ 0.1 -4 0.3 at the scale of rich clusters, and its evolution with time 
parametrized in terms of the parameters rj ~ 4- 0.9 and s ~ 4- 1.5) are 
discussed to point out the role of different, possible gas-phase phenomena in 
cluster evolution. Here we take into account the minimal and maximal levels of 
IGM evolution that are still consistent - given the observational uncertainties 
- with the current data ([65] [25] [26] [28]). 

So, in Fig. 8 we show the possible range of J 7 produced by a population of 
galaxy clusters considering reasonable uncertainties in the theoretical descrip- 
tion of the clusters structure and of their distribution in mass for the consid- 
ered cosmological models. According to our previous estimation of the uncer- 
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tainties, it is reasonable to expect that galaxy clusters could provide up to a 
few % of the EGRET [1] DGRB. Higher values of I y would require a choice 
Lcr ^ 10 44 erg/s, or a large contribution from galaxy groups: both these 
possibilities seem to be quite unrealistic. 



7 High energy neutrinos from clusters of galaxies 

Gamma rays from galaxy clusters are always accompanied by neutrino produc- 
tion due to the decay of charged pions, according with eq.(2). In Fig. 9 we show 
the range of cluster contribution to the diffuse neutrino background (DNB), 
I u , as a function of the neutrino energy E, expected in the viable models for 
structure formation. We considered in the predictions of Fig. 9 the level of the- 
oretical uncertainties in the description of the cluster structure and evolution 
(see Sect. 6 for a discussion). We considered the whole contributions due to the 
fluxes of and u e , and we compare them to the level of the atmospheric neu- 
trinos as given by Gondolo et al. [77] . The diffuse v flux produced by clusters 
is at a level of F u ^ 1CT 2 neutrinos km~ 2 yr^sr^GeV^ 1 at E ^ 3- 10 5 GeV, 
where the spectral shape of the cluster fluxes (oc E~ 2A ^ 2 - 3 ) emerges from the 
E distribution of the atmospheric z/'s (see Fig.9). However, such an estimate 
has to be considered quite optimistic because it refers to the specific choice 
of parameters that maximize J 7 . In such a maximal model for 7 7 , at E — 10 3 
GeV one could observe a sensitive diffuse flux from galaxy clusters at a level 
J 7 ~ 30 4- 2 • 10 3 neutrinos km~ 2 yr~ A sr^GeV' 1 . This diffuse emission could 
be one of the most intense neutrino backgrounds of cosmological origin at 
these energies, comparable only to that produced by unresolved AGN cores 
[71]. Single clusters could be detected with the next generation, large neu- 
trino telescopes as the signals from rich, nearby Coma-like clusters are of the 
order F v m 10 3 4 10 4 neutrinos kmr 2 yr~ l (at z ^ 0.05) at E ~ 10 3 GeV. 
Much distant clusters are more difficult to be detected as their fluxes scale as 
oc d L 2 (z,tt ) [see eqs.(34-35) and Fig.2]. 



8 Discussion and conclusions 

In this paper we presented a detailed study of the diffuse emission of 7-rays 
and neutrinos from clusters of galaxies. Using realistic modelling of the cluster 
structure, of their formation history and of their evolution with cosmic time, 
we found that galaxy clusters can provide ^ 1% of the DGRB measured by 
EGRET (in the first release by OWZ [1]). 

Our estimate of J 7 is quite independent on the geometry of the universe, 
on the assumed cosmological model and on the amount of IGM evolution, 
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because most of their contribution to the DGRB comes from nearby, z ^ 0.2, 
clusters. In fact, at these redshifts the effects of curvature do not take place 
strongly in changing the perturbation growth factor, V(z,Q ), (normalized 
at the present epoch), the difference in cluster evolution are small when the 
different models are normalized to the local abundance of clusters observed 
in X-rays and the available amount of IGM evolution - even if considered to 
be quite strong, f g oc (1 + z) _1 ^~ 2 - can provide only small variations to the 
cluster 7-ray luminosities, as L 7 oc f g (see eq.18). On account of all these 
aspects, we consider that our results for the contribution of galaxy clusters to 
the DGRB are quite robust. 

Our approach differs substantially from the previous ones in several (among 
others) aspects: 

i) we considered - differently from all the previous approaches - a self-consistent 
approach to the formation of clusters following the spherical collapse model 
[69] complemented with a realistic IGM density profile, consistent with the 
most recent determinations from X-ray observations. This fact has important 
effects on the CR confinement within cluster cores and hence on the relative 
7-ray and neutrino emission rates; 

ii) we considered (as BBP did) here an energy dependent diffusion coefficient 
which results in a very general picture of the CR confinement within cluster 
cores; 

in) we also considered here - at variance with the previous approaches - the 
effects of a possible evolution in the cluster IGM content. This is consistent 
with the present indications of a variation in the IGM content from groups to 
rich clusters in the local frame and with the X-ray, shock (or entropy) induced, 
luminosity evolution observed from numerical simulations [42] and predicted 
in analytical models (both shock and entropy models) for the evolution of 
X-ray clusters [48] [40]; 

iv) we use the PS cluster MF that was found to be consistent with N-body sim- 
ulations over a large dynamical range and up to z ^ 2 [72] . We normalized it to 
the local abundance of clusters detected in X-rays. In the previous approaches 
average values for the overall cluster abundance, n c \ ~ 4 4- 7 • 10~ 5 Mpc~ 3 , 
were used [11] [13] without considering the effect of an evolving cluster mass 
function; 

v) using a self-consistent modelling of the IGM we found an analytical cor- 
relation between 7-ray and X-ray emission for clusters. The predicted ratio 
F^/Fx oc f g 1 r~ 1 T~ 1 / 2 [see eq.(32)] provides a behaviour of the F 7 — Fx rela- 
tion different from that obtained by Ensslin et al. [20], F 1 /Fx oc T 1//2 , because 
we did not assume any (partial) equipartition between IGM thermal energy 
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and relativistic jet particles. Our correlation results only from the basic elec- 
tromagnetic and hadronic emission mechanisms in which the IGM protons are 
the targets for both the X-ray thermal bremsstrahlung emission and for the 
pp collisions responsible for 7-rays. 

vi) using such a correlation we derived a sample of nearby clusters with pre- 
dicted 7-ray fluxes observable with the next generation 7-ray telescopes. Inci- 
dentally, we found a 7-ray flux for Coma, F 7 (> lOOMeV) ~ 8.5 ■ 1CT 9 photons 
s _1 cm~ 2 which is consistent with the EGRET upper limits for this clus- 
ter (previous specific predictions [12] [20] seem to exceed the EGRET upper 
limit). 

Our numerical results for 1 1V are in reasonable agreement with those obtained 
by BBP, even though based on a quite different description of the cluster 
structure and evolution. This agreement is due to the fact that BBP considered 
a constant comoving density of clusters, n c \ ~ 5 • 10~ 5 Mpc~ 3 , assumed to 
be a fair sample of the baryons in the universe, and containing a fraction 
Q b ~ 0.5Qbbn (where flsBN is the value of the baryon density predicted by 
Big Bang Nucleosynthesis). Under these assumptions, BBP obtained a value 
for J 7 higher by a factor ~ 3 -r 4 with respect to our result, based on values 
f g ~ 0.1 (see Sect. 2). Our refined calculations show why their assumption of 
considering n c i ~ const was reasonable: the clusters effectively contributing 
to the DGRB are located at z ^ 0.2 (see Fig.5), where the effects of evolution 
do not have room to take definitely place (see Fig. 7). 

Because of the inherent uncertainties in the predictions of quantities whose cal- 
culation involve to set the values of parameters which are not known precisely, 
we also estimated the range spanned by I 1)U for the combination of parame- 
ters allowed by the present observational ranges. In fact, the description of the 
cluster structure and evolution that we used in our analytic approach consider 
only ensemble averaged quantities. But we observe a whole distribution of the 
real cluster properties with respect to the average cluster moulding. Some 
amount of variance is needed to be considered in cluster modelling in order to 
ensure the predictive power of the viable models for structure formation. To 
explore the role of the uncertainties in the relevant quantities we considered 
several sources of uncertainties (see Sect. 7). 

From an inspection of Figs. 9 and 10 we note that the effects of the possible 
theoretical uncertainties in the description of the cluster and IGM evolution 
could change the predicted contributions for the DGRB and for the DNB by 
a factor ^ 3, setting the maximal level of 7 7 to a few % of the EGRET value. 

The DGRB seems to be mostly produced by AGNs (FSRQ and/or BLLacs) 
and/or blazars (we take here an estimate of ~ 60 -7- 65% [73] of the EGRET 
diffuse flux [1]). Diffuse 7-ray emission could be observed also from normal 
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galaxies yielding a contribution ~ 5% [73]. When added to the ~ 10 4- 15% of 
the DGRB contributed by their high-i? photons interacting with other existing 
backgrounds (e.g. IR, CMB [73]) one gets only ~ 15 4 20% of the DGRB left 
for truly diffuse or extended sources. Of this amount, a fraction of the diffuse 
7-ray flux ~ 3 4- 5% is predicted [74] to originate from decaying topological 
defects (see [75]) and interactions of UHE particles with the CMB. Note, how- 
ever, that the amount and the spectral distribution of this possible diffuse 
background depend sensitively on the amplitude of the primordial magnetic 
field on scales larger than supercluster sizes. So, according to the previous 
estimates, the presence of all these sources of diffuse 7-ray emission (even 
though partially model dependent) determines an upper limit to the contri- 
bution of extended extragalactic sources to the DGRB, that is ~ 10 4- 22% of 
the OWZ EGRET level [1]. This sets rather weak constraints on the level of 
CR production in clusters and hence on the presence and activity of AGNs in 
clusters or on the formation and efficiency of accretion shocks around clusters. 
However, if we consider the revised level of the DGRB as derived by SWZ [2], 
then the previous upper limit reduces to ^ 2%. Our predictions of the DGRB 
contributed by galaxy clusters I 1;C i/^egret ~ 0.005 4 0.02 is perfectly com- 
patible with the presence of both a population of evolving FSRQ and AGNs 
dominating the 7-ray sky and with the presence of truly diffuse backgrounds 
like those previously discussed. Note, however, that the major source of uncer- 
tainty in the level of the extragalactic DGRB comes from the contribution of 
the AGNs. A fluctuation analysis of the EGRET data is needed to have more 
definite indications on the level of the DGRB contributed by discrete sources. 
If, on the other hand, CR acceleration will be found to be relevant in clusters 
(yielding Lqr substantially larger than 10 44 erg/s), then the predicted level 
of DGRB produced by galaxy clusters can set interesting limits to the space 
density and evolution of 7-ray AGNs. 

The sensitivities and angular resolutions achievable by the next generation 
gamma-ray (INTEGRAL, GLAST, AMS) and neutrino (see [76] for a review) 
detectors will be able to shed a new light on the high energy phenomena 
occurring in large scale structures. 
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Table 1. A list of predicted 7-ray fluxes from clusters of galaxies. 
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Table caption 

( a ) n c in units 10~ 3 cm -3 . 

( 6 ) X-ray fluxes in the (2 — 10) keV band in units 10~ 10 erg cm" 2 s _1 . 
( c ) Gamma ray fluxes F 7 (> lOOMeV) in units 10~ 8 photons s" 1 cm" 2 
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Fig. 1. Gamma ray luminosities for clusters in different cosmological scenarios: 
CDM (Qq = 1; dashed curve), CDM+A (Uq = 0.4; continuous curve), open CDM 
(£Iq = 0.3; long-dashes curve) and mixed DM (Qq = 1, VL U = 0.3 ; dotted curve). 
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Fig. 2. Gamma ray fluxes for clusters at z = 0.023, the redshift of the Coma cluster. 
The arrow indicates the EGRET upper limit for A1656 (Coma). Curves labelled as 
in Fig.l. 
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Fig. 3. Cluster 7 -ray fluxes for the objects listed in Table 1 plotted against their IGM 
temperature (filled dots). The curves represent the F^ — T relationship in a CDM+A 
(Qo = 0.4, h = 0.6) cosmology at redshift z = 0.023 (continuous curve) and at the 
maximum redshift of the clusters in our sample (dashed curve). 
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Fig. 4. The frequency of clusters with a predicted j-ray flux in the sample of Table 
1. 
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Fig. 5. The redshift distribution of clusters with fluxes F^(> WOMeV) in a CDM 
model with (rj = 0.2, s = 1.8: continuous histogram) and without (rj = 0.2, s = 0: 
dotted histogram) IGM evolution (see text for details). 
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Fig. 6. The contribution of galaxy clusters to the DGRB in different cosmological 
models: CDM (Qq = 1; dashed curve), CDM+K (Uq = 0.4; continuous curve), open 
CDM (Qq = 0.3; long-dashes curve) and mixed DM (Q,q = 1, Q, v = 0.3 ; dotted 
curve). 



34 



1 I I III 




0.1 1 10 100 

Log E (GeV) 



7. The predicted level of DGRB from galaxy clusters in a CDM+A 
= 0.4, h = 0.6) model with (rj = 0.2, s = 2.2: continuous curve) and without 
0, s = dashed curve) IGM evolution (see text for details). 
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Fig. 8. The expected DGRB from clusters, considering various sources of theoretical 
uncertainties in the cluster modelling. A flat CDM (Qo = l;h = 0.5; n = I) model 
is considered here. 
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Fig. 9. Same as Fig. 8 but for the diffuse v background. Heavy dashed line show the 
diffuse flux from atmospheric v's (Gondolo et al. 1995). 
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